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ABSTRACT 


A  method  for  the  calculation  of  turbulent  shear  flows  by 
an  invariant  modeling  closure  of  the  equations  for  the  second- 
order  correlations  of  fluctuating  quantities  in  a  turbulent 
medium  is  described.  The  relationship  of  the  primary  computa¬ 
tional  scale  used  in  this  model  and  the  longitudinal  integral 
scale  of  the  turbulence  under  consideration  is  determined.  With 
these  results  in  hand,  a  technique  is  described  by  which  it 
should  be  possible  to  determine  the  intensity  and  integral  scale 
of  turbulence  measured  by  instrumented  aircraft  when  only  partial 
spectra  are  available.  This  is  possible  when  the  distributions 
of  mean  wind  velocity  and  temperature  in  the  atmosphere  in  the 
vicinity  of  the  tests  are  available.  Preliminary  calculations 
using  the  method  for  a  case  when  both  the  intensity  and  integral 
scale  of  turbulence  were  known  yielded  very  good  agreement 
between  computed  and  experimental  results. 
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INTRODUCTION 


x  . 


Some  time  ago,  the  senior  author  of  this  paper  proposed  an 
approach  to  the  computation  of  turbulent  shear  flows  (Ref.  1) 
which,  because  of  it;  generality,  gave  promise  of  allowing  one 
to  make  calculations  of  turbulent  atmospheric  and  vortex  motions. 
In  1969  this  method  in  its  earliest  form,  was  applied  to  the 
calculation  of  the  generation  of  turbulence  in  the  atmosphere 
(Ref.  2).  The  results  of  these  calculations  were  intriguing 
enough  to  warrant  refinement  of  the  method  and  its  further  appli¬ 
cation  to  calculations  of  the  generation  of  atmospheric  turbulenc 
During  fiscal  year  1971  support  was  obtained  from  two  sources 
towards  these  two  ends.  A  moderate  level  of  support  was  obtained 
from  NASA  under  Contract  No.  NASW-1868  for  further  refinement  of 
the  method  and  a  somewhat  higher  level  of  support  was  obtained 
from  t lie  Air  Force  for  development  cf  the  method  for  application 
to  a  specific  problem  faced  by  the  Flight  Dynamics  Laboratory; 
this  report  discusses  the  results  obtained  for  tie  Air  Force. 


II.  STATEMENT  OP  THE  PROBLEM 


A  problem  that  has  faced  and  continues  to  face  the  Flight 
Dynamics  Laboratory’s  Design  Criteria  Branch  is  that;  of  deter¬ 
mining  two  fundamental  characteristics  of  atmospheric  turbulence 
under  any  given  set  of  conditions.  These  characteristics  or 
parameters,  which  define  in  a  gross  way  the  structural  response 
of  aircraft  to  atmospheric  turbulence,  are  the  root  mean  square 
of  the  vertical  component  of  the  turbulent  velocity  field  o„ 
and  an  integral  scale  of  the  turbulent  velocity  field  L  .  In 
general,  it  is  assumed,  for  work  concerned  with  structural  problems 
that  the  turbulent  velocity  field  is  sufficiently  isotropic  so 
that  only  two  integral  scales  must  be  considered.  One  is  the 
longitudinal  integral  scale  L  defined  by 

syco 

L  =  L, ,  =  - - -  I  <u’(x)u.(x  +  x,  )>dx,  (1) 

11  <u’u’Xx)Jo  1  “  1  “  11 

The  other  is  the  transverse  integral  scale  L^ 

/?« 

L_o  =  - - -  I  <u’(x)u’(x  +  xn)>dx  =  L?9  (2) 

33  <u^u^>(x)  J o  3  3  1  1  22 

It  is  well-known  from  the  theory  of  isotropic  turbulence  (Ref.  3) 
that  these  scales  are  related;  namely. 


L  =  L11l  =  2L33  =  2L22  (3) 

For  the  calculations  involved  in  practical  structural  design,  it 
has  been  assumed  that  Eq.  (3)  is  valid,  so  that  the  problem  facing 
the  structural  designer  is  that  of  determining  a  and  L  for 
any  given  atmospheric  situation.  * 

In  principle,  the  determination  of  these  two  parameters  is 
straightforward.  In  practice,  however,  the  task  is  impossible  in 
many  cases  because  accurate  measurements  of  turbulence  intensity 
are  not  available  over  a  wide  enough  range  of  wave  numbers. 
Consider  the  spectrum  of  turbulence  shown  in  Figure  1.  This 
spectrum,  taken  from  Ref.  4,  represents  the  power  spectral 
density  <1 >w(ft)  of  the  turbulent  vertical  velocity  w’  measured 
in  a  thunderstorm,  plotted  against  the  reduced  frequency  Q  . 

This  spectrum  is  typical  of  many  atmospheric  spectra  that  have 
been  obtained.  The  trustworthy  portion  of  the  spectrum  exhibits 

a  <t>  —  n"5/3  behavior  throughout.  This  behavior,  while  it  is 

gratifying  to  students  of  turbulent  motions,  does  not  permit  one 


2 


. n.  ’  N~nT'  < -^r^.’;^«rr*^»r,r^’iT»^r^rr 


to  estimate  either  the  integral  scale  or  the  intensity  of  the 
turbulence.  It  has  been  pointed  out  by  Houboit  that,  if  the 
actual  spectrum  i_  a  von  karman  spectrum,  namely, 

r£h  1  +  4  (1.539 La)2 

T,.  =  -4 - * - tt-tt-7-  (*0 


w 


V 


LI  <■  [l.V'^L 


] 


one  can  infer  from  a  record  such  as  that  shown  in  Figure  1  only 
the  following  relationship  between  the  intensity  and  the  integral 
scale : 


o2  -  -  1 
w  -  0.521 


W 


(5) 


where  fi*  ,4>$  represents  any  point  on  the  minus  five-thirds  slope 
portion  of  the  spectrum.  Two  spectra  satisfying  Eq .  (A)  are  shown 
in  Figure  1.  One  spectrum  is  for  o w  =  32.33  ft/sec  and  L  =  5600 
feet,  while  the  second  is  for  ow  =  25.66  ft/sec  and  L  =  2800  feet. 

This  example  illustrates  the  difficulty  of  finding  the 
quantities  ow  and  L  when  only  partial  spectra,  such  as  the  one 
just  discussed,  are  given.  The  study  supported  under  the  present 
contract  was  an  effort  to  develop  a  method  by  which  it  might  be 
possible  to  estimate  both  ow  and  L  from  partial  spectra 
when  the  profiles  of  mean  wind  and  atmospheric  temperature  were 
available . 

The  technique  proposed  was  the  following.  The  method  of 
invariant  modeling  permits,  as  is  shown  below,  the  computation  of 
the  turbulent  structure  of  the  atmosphere  if  the  local  mean  wind 
and  temperature  profiles  are  known.  This  computation  is  possible 
if  a  scale  parameter  Aq  ,  related  to  L  ,  is  assumed.  The 
relationship  between  computed  0  and  assumed  L  is  different 
than  that  given  in  Eq.  (5).  Thul ,  if  the  relationship  of  L  and 
Aq  is  known,  it  is  possible  to  determine  both  and  ow  by 

finding  the  intersection  of  the  curve  representing  Eq.  (5)  and 
of  the  curve  c btained  from  the  results  of  the  computations. 

The  study  undertaken  for  the  Air  Force  consisted  of  two 
parts.  First,  detailed  investigations  were  carried  out  in  which 
computations  of  turbulence  in  a  free  shear  layer  and  in  an  axially 
symmetric  free  jet  were  made.  In  these  studies,  the  local  value 
of  A^  was  found  In  each  case  that  best  reproduced,  within  the 

ability  of  the  mathematical  model,  the  turbulent  structure  of  the 
particular  flow  under  investigation.  This  was  possible  because 
of  the  fairly  large  amount  of  detailed  information  concerning  the 
turbulent  structure  of  these  two  flows.  The  local  value  of  Aq 


(ft/sec)  /(radians /ft) 


£2,  radians/ft 

Fig.  I.  Power  spectrum  measured  in  a  thunderstorm 
(Ref.  4)  and  two  von  Kdrman  spectra  fitted  to  it 


was  then  compared  with  the  local  value  of  L  as  measured  in  each 
case.  In  this  way,  it  was  determined  that  the  following  relation¬ 
ship  between  L  and  is  representative: 

Ax  -  0 . 6L 

Second,  with  this  result  in  hand,  a  computation  was  carried  out, 
as  an  initial  test  case,  to  see  if  the  method  discussed  above  for 
recovering  atmospheric  scales  and  intensities  was  plausible. 
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III.  DESCRIPTION  OP  AN  INVARIANT  MODEL  OP  TURBULENT  SHEAR  LAYERS 


The  equation  for  the  time-independent  mean  velocity  in  an 
incompressible  turbulent  medium  was  given  many  years  ago  by 
Reynolds  (Ref.  5).  It  is 


puJu 


,5i,j  =  'Li  +(^  -  p<uJ'uP),j 


In  this  paper,  bars  over  a  quantity  or  angular  brackets  around 
a  quantity  indicate  average  values  of  that  quantity  while  primes 
indicate  the  instantaneous  fluctuation  of. the  quantity  from  its 
mean  value.  The  mean  molecular  stress  is  given  by 

* l  -  gJk|1(3i,k +  \,i)  (7) 

The  second-order  correlation  of  the  velocity  fluctuations  that 
appears  in  Eq.  (6)  represents  the  transport  of  momentum  by  turbulent 
eddies  and  is  called  the  Reynolds  stress.  An  equation  for  this 
second-order  tensor  was  also  given  by  Reynolds.  It  is,  for  steady 
flow , 

p"J<uiu>>,j  -  -p<uJ'uiP"i,j  -  p<uJ'ul>"k,j  -  p<uJ'uiui?,j 


-  <uip'\k  -  <ukP’>,t  *  <p‘(ui,k  *  «i,i» 


+  ;'.g  ‘7uijmukfD 


In  the  past,  it  has  been  customary  to  carry  out  investigations 
of  turbulent  shear  flows  by  means  of  Eq.  (6).  In  these  studies 
the  unknown  second-order  velocity  correlation  term  was  modeled 
in  terms  of  the  mean  velocity  and  a  length.  Computations  made  in 
this  manner  form  the  vast  bul.c  of  the  literature  on  turbulent  shear 
flow  calculations  tc  the  present  time.  The  methods  that  are  now 
in  use  employing  this  type  of  modeling,  having  evolved  over  a  period 
of  many  years,  are  exceedingly  useful  and  enough  empirical  data 
have  been  amassed  to  enable  ono  to  calculate  solutions  to  a  wide 
variety  of  engineering  problems  with  a  great  deal  of  confidence. 
There  are,  however,  a  number  of  problems  facing  engineers  today 
which  require  the  calculation  of  turbulent  shear  flows  for  which 
there  is  nc  existing  body  of  experimental  data.  Two  flows  which 
come  readily  to  mind  are  the  generation  of  turbulence  and  turbulent 
transport  correlations  by  the  earth's  atmosphere  and  the  decay 
of  a  turbulent  vortex.  In  the  case  of  these  two  flows,  we  may  ask 
ourselves  the  following  questions:  "Is  there  not  a  somewhat  more 
basic  method  of  computing  turbulent  transport  phenomena  than  the 
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eddy  viscosity  or  mixing  length  models  presently  in  use?"  "Can 
not  such  a  method  permit  us  to  generalize  the  experimental  results 
that  presently  exist  so  as  to  estimate  the  nature  of  turbulent 
flows  that  have  not  yet  been  investigated  experimentally?"  The 
answers  to  these  two  questions  are  not  as  straightforward  as  one 
would  like.  In  answer  to  the  question  as  to  whether  there  is  a 
more  basic  method,  the  reply  must  be:  Yes,  but  the  real  difficulty 
comes  when  one  tries  to  establish  just  how  much  more  fundamental 
the  proposed  new  method  is  to  be.  If  a  new  method  is  truly  more 
fundamental,  then  it  will  allow  better  estimates  of  unknown  flows 
than  can  be  made  by  older  techniques.  It  is  fair  to  say,  at  the 
present  time,  that  there  is  considerable  hope  among  those  who  practice 
the  art  of  calculating  turbulent  shear  flows  that  the  new  methods 
now  being  developed,  of  which  the  method  to  be  described  here  is 
but  one,  take  into  account  enough  physics  that  is  not  contained 
in  older  models  so  that  a  technological  advance  can  be  made.  It 
is  as  yet  too  early  in  the  history  of  these  new  methods  to  make  any 
strong  statement  as  to  just  how  much  more  powerful  they  are  than 
the  older  methods. 

If  one  wishes  to  make  use  of  both  Eqs .  (6)  and  (8)  in  computing 
turbulent  shear  flows,  the  first  step  must  be  a  choice  of  models 
for  those  terms  in  Eq.  (8)  which  are  not  expressed  in  terms  of  the 
mean  velocity  or  the  second-order  velocity  correlation.  The  terms 
which  must  be  modeled  are 

(1)  the  pressure-strain  correlation  in  the  tendency-towards- 
isotropy  term,  namely,  <p'(u£  k  +  u£  i)>  ; 

(2)  the  third-order  tensor  in  the  velocity  diffusion  term, 
namely,  <ujuju£>  ; 

(3)  the  pressure  velocity  correlation  in  the  pressure  diffusion 
terms,  namely,  <(  ujp  O  ; 

(4)  the  general  viscous  dissipation  term  2|xgmn<(u'  u,'  3 

1  y  m  K  jil 

There  are  many  ways  in  which  a  modeling  of  the  above-noted 
terms  may  be  accomplished.  We  have  tried,  for  our  initial  investiga¬ 
tion,  to  take  as  simple  a  model  as  possible  for  each  term.  We 
have  then  attempted  to  determine  by  calculation  the  adequacy  of 
the  chosen  model  and  the  sensitivity  of  the  calculated  results 
to  the  particular  choice  of  model. 

To  date  the  following  models  have  been  investigated  to  some 
extent . 

(1)  For  the  tendency-towards-isotropy  term,  we  choose,  following 
Rotta  (Ref.  6),  the  following  model: 


<»’<uU  *  ut,i)>  ’  '  ^(<uiufc>  -  snc§) 


(9) 
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i 


where 


q 


2  „  K  =  <um'u,p  =  <u'2>  *  <v'2>  *  <»'‘> 


I2" 


(10) 


and  A]_  is  a  scalar  length  associated  with  the  tendency  towards 
isotropy  and  is  to  be  identified.  More  complicated  models  of 
this  term  have  been  discussed  by  Chou  (Ref.  7)  and  by  Hanjalic 
and  Launder  (Ref.  8).  To  date,  we  have  had  considerable  success 
using  the  simple  Rotta  model  given  above  and,  in  line  with  our 
philosophy  of  using  the  simplest  possible  model  that  will  give 
reliable  results,  we  have  confined  the  majority  of  our  computa¬ 
tional  studies  to  the  use  of  Eq  (9). 


(2)  For  the  velocity  diffusion  term,  we  must  model  <(u,'u '.ul  )>  . 
The  simplest  covariant  tensor  of  rank  three  that  is  symmetric  inK 
all  three  indices  that  we  can  form  out  of  the  second-order 
correlations  's 


<uiuj>,k  *  <ujuk>  i  +  <uiul>  j 


This  expression  has  all  the  tensor  and  symmetry  characteristics 
required  of  our  model.  To  make  it  dimensionally  correct,  the  above 
expression  must  be  multiplied  by  a  scalar  velocity  and  a  scalar 
length.  The  simplest  scalar  velocity  we  can  form  from  the  second- 

order  correlations  is  /  <(  um  u’N  =  q  ,  so  we  model  the  tensor 

.  .  N  in  ^ 


<  UiUjUk  ^  as 


<uiujui>  ■  -V  <uiujlk  +<ujuk>  i  -<UK >,} 


(n) 


where  Ay  is  a  scalar  length  associated  with  velocity  diffusion 
and  is  also  to  be  identified  by  matching  experimental  results. 


(3)  The  pressure  velocity  correlation  <(p'u^)>  in  the  pressure 


diffusion  term  is  modeled  by  analogy  with  the  velocity  diffusion 
term  as 


<P'U;£>  “  -pA3q<ura'u^> 


m 


(12) 


In  our  work,  in  order  to  cut  down  the  number  of  parameters  in  our 
turbulence  model,  we  have  considered  only  two  special  cases  of 
Eq .  (12).  We  have  considered  the  case  A,  =  A2  and  the  case  where 
<p'Ujp>  =  0  ,  i.e.,  the  case  A^  =  0  .  ^ 


(i|)  We  have  considered  two  models  for  the  expression 


g  n  <C  u!  ui  0  appearing  in  the  viscous  dissipation  term: 

x  j  til  k  j  n 

<uiuk> 


(a) 


gmn<u 


i,m  k,n 


>- 


(13) 
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(14) 


«ysj?,F»r»T!SEw  ■ 


and 


(b) 


iTCuj 

X  > 


muk,n'> 

/ 


In  both  these  models,  Hs  a  dissipative  length  scale.  The  argu¬ 
ment  for  choosing  the  latter  expression  is  that  it  is  expected 
that  the  turbulence  will  be  almost  isotropic  in  that  part  of  the 
spectrum  responsible  for  dissipation  of  turbulent  kinetic  energy. 
Thus,  one  would  expect  the  dissipation  to  be  almost  isotropic  even 
if  the  turbulence  itself  is  not.  Further,  there  is  experimental 
evidence  that  the  loss  of  shear  correlation  by  viscous  action  is 
relatively  much  smaller  than  the  loss  of  kinetic  energy  by  viscous 
action.  In  our  initial  computations  using  Eq.  (1*1)  as  a  model  of 
dissipation,  we  experienced  some  difficulties  in  obtaining  solu¬ 
tions.  There  was  a  tendency  for  solutions  to  develop  with  negative 
values  for  the  mean  square  velocities  when  the  turbulence  became 
very  nonisotropic.  This  tendency  was  overcome  by  the  use  of  Eq . 
(13)  for  the  dissipation  model.  This  model  does  not  have  a  large 
effect  on  the  development  of  the  shear  correlations  because  the 
primary  contribution  to  loss  of  shear  with  this  model  is  not  the 
dissipation  term  but  the  tendency-towards-isotropy  term.  Although 
the  whole  question  of  modeling  the  dissipation  term  is  still  under 
investigation,  the  work  reported  here  was  carried  out  using  Eq. 
(13),  for  the  reasons  stated  above. 

In  the  models  given  above,  we  would  expect  that  the  scalar 
lengths  A^,  and  would  all  be  related  to  the  local  integral 

scale  of  the  turbulence.  These  lengths  are,  in  turn,  related  to 
the  local  scale  of  the  mean  motion  for  the  flows  we  shall  investi¬ 
gate  here,  and  we  make  the  assumption  in  the  computations  we  will 
discuss  presently  that  A^,  A^,  and  A^  are  all  proportional 

to  some  local  characteristic  length  6  °f  the  mean  motion 

under  consideration. 

We  will  expect  the  length  appearing  in  the  dissipation  model 
to  be  related  to  the  microscale  of  the  turbulence  which,  in  turn, 
must  be  related  to  the  integral  scale  via  a  Reynolds  number  in 
such  a  way  that  production  of  turbulence  is  balanced  to  a  large 
extent  by  dissipation. 

If  the  models  we  have  just  discussed  [Eqs.  (9)  and  (11)~(13)] 
are  substituted  in  the  basic  equation  for  the  second-order  velocity 
correlations  [Eq.  (8)],  the  resulting  equation,  taken  together 
with  the  momentum  and  continuity  equations,  makes  a  closed  set 
(see  Appendix  A).  When  this  set  is  reduced  to  boundary  layer  form, 
it  is  found  to  form  a  parabolic  system  (see  Appendix  R). 
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I'Tv” 


This  set  of  equations  will  admit  similarity  solutions  at 
high  Reynolds  numbers  as  well  as  permit  calculations  of  turbulent 
flows  near  walls,  if  one  makes  the  following  choice  of  the  relation 
between  the  length  scales: 


and 

where 


> 

H-1 

II 

clAchar 

(15) 

II 

OJ 

< 

G2A1  =  c26char 

(16) 

A3  _ 

c3A1  *  °36ehar 

(IT) 

X  = 

A./va  +  b  •  Re. 

1  Ai 

(18) 

ReA,  = 

pqA^p. 

(19) 

For  self-similar  free  turbulent  flows,  the  structure  given 
above  is  all  that  is  needed  to  compute  a  turbulent  shear  layer 
or  a  free  jet,  provided  the  five  constants,  c^ ,  c2,  c^,  a  ,  and  b 

are  given.  To  find  these  constants,  we  must  resort  to  the  compari¬ 
son  of  calculated  flow  fields  with  experimental  results. 


If  we  wish  to  compute  a  boundary  layer  flow,  we  must  consider 
an  additional  problem.  When  a  wall  is  present  in  a  shear  flow, 
we  wish  to  apply  the  boundary  condition  at  the  wall  that 


<Cuiu,'^>  =  0 

xiY  z=o 

where  z  is  measured  normal  to  the  surface.  In  addition,  there 
should  be  no  diffusion  of  through  the  surface,  so  that 

S  <(  u^Ujp’/dz  -  0  at  z  =  0  .  Thus,  it  is  reasonable  to  assume 

that  near  the  wall 

■  Atkzltn  <20> 

where  A ^  is  a  constant  and  Tj  is  a  positive  constant.  But  if 

there  is  no  diffusion  through  the  well,  then  all  that  is  diffused 
towards  the  wall  by  viscosity  at  z  =  e  is  dissipated  in  the  region 
between  z  *  €  and  z  =  0  ,  (It, is  easily  verified  that  all  other 
terms  in  the  model  equation  for  <!u!ul/  are  neglible  if  e  is 
small  enough.)  Thus,  1  K 


<uiuk> 


dz  =  u 
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or,  using  Eq.  (20) 


(1  +  T))er‘ 


Thus,  near  a  solid  surface,  we  will  always  assume,  in  applying  our 
model,  that  Eq.  (21)  holds  in  the  region  near  the  wall. 

It  is  convenient  to  express  this  result  in  terms  of  A-,  . 

Near  a  wall,  Eq.  (18)  becomes 

X  =  A1/v/a  (22) 

Using  Eq.  (21),  we  may  write 

A-^  =  a\/a  z  (23) 

Thus,  for  boundary  layer  flows,  a  is  another  number  which  must 
be  found  from  experimental  results. 

In  our  first  attempts  to  construct  a  model  of  turbulent  shear 
flows  (Refs.  9  and  10),  the  following  assumptions  were  made  to 
construct  the  simplest  possible  model  of  boundary  layer  flows: 

(1)  It  was  assumed  that  all  the  large  lambdas  associated 
with  inviscid  modeling  were  equal,  i.e.,  A^  =  Ag  =  A^  =  A  . 

(2)  It  was  assumed  that  a  was  equal  to  one. 

(3)  In  the  outer  portion  of  a  boundary  layer,  A  was  taken 

to  be  a  constant  c-^  times  <5  (6  is  the  value  of  z  for 

which  u  is  99%  of  the  free  stream  velocity).  This  value  was 
assumed  to  hold,  independent  of  z  ,  as  the  wall  was  approached, 
until  A  became  equal  to  /a  times  z  .  For  smaller  values  of  z  , 
A  was  taken  equal  to/a  z  , 

With  these  assumptions,  the  boundary  layer  forms  of  Eqs.  (6) 
and  (8)  with  appropriate  modeling  (Appendix  B)  were  solved  with 
various  choices  for  the  parameters  a  ,  b  ,  and  c^  =  A/6  ^  to 

produce  a  developing  turbulent  boundary  layer  on  a  flat  plate. 
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It  was  determined  at  that  time  that  the  following  choice  of 
parameters 


C1  =  A/6.99  =  0,064 

a  =  2.5  (24) 

b  =  0.125 

yielded  a  fair  representation  of  a  turbulent  boundary  layer.  The 
mean  velocity  profile  and  the  behavior  of  skin  friction  with 
Reynolds  number  were  adequately  represented.  The  distributions 
of  the  second-order  correlations  within  the  boundary  layer  were 
reasonable . 

The  results  of  this  original  parameter  search  were  used  to 
compute  a  number  of  other  turbulent  flows  in  order  to  demonstrate 
the  method  (Refs.  9  and  10). 

Before  proceeding  with  further  applications,  it  was  considered 
necessary  that  a  more  detailed  parameter  search  should  be  made. 

In  particular,  two  free  turbulent  flows  -  the  free  jet  and  the 
free  shear  layer  -  should  be  calculated  to  determine  the  values  of 
the  parameters  c?  ,  c,  ,  a  ,  and  b  that  would  best  fit  the  exper¬ 
imental  results  for  both  flows.  (The  equations  for  the  free  shear 
layer  are  the  same  as  those  for  the  boundary  layer  given  in  Appendix 
B.  The  equations  for  the  axially  symmetric  free  jet  are  given 
in  Appendix  C.)  The  value  of  c,  being  the  ratio  of  A,  to  some 
arbitrarily  defined  characteristic  length  in  each  case  is  not  an 
invariant  of  the  problem  and  was  to  be  chosen,  with  fixed  values 
of  the  other  parameters,  to  obtain  best  results  in  each  ease. 

Once  these  studies  were  complete,  the  model  would  be  used  to  compute 
turbulent  boundary  layer  flows  so  that,  by  comparison  with  experimental 
results,  values  for  c,  and  a  could  be  made  for  this  flow. 

Hopefully,  all  flows  could  be  described  in  a  reasonable  way  by  a 
single  choice  of  the  basic  model  parameters  c-  >  ,  a  ,  b  , 

and  (where  appropriate)  cx.  The  values  of  local  A^-5  determined 
from  the  values  of  c,  in  each  case  were  then  to  De  compared 
with  the  local  magnitude  of  the  integral  scale  l  in  each  case. 

If  it  was  found  that  the  value  of  c,  represented  a  choice  that 
amounted  to 

A  ®  const  L  ®  PL  (25) 

then  it  would  be  assumed  that  a  reasonably  invariant  model  had 
been  determined. 


IV.  THE  SEARCH  FOR  NEW  MODEL  PARAMETERS 


Our  search  for  a  new  model  of  turbulent  shear  layers  began 
with  an  attempt  to  describe  the  axially  symmetric  free  jet  with 
the  original  turbulence  model  obtained  for  a  boundary  layer  flow. 
This  model,  as  mentioned  in  the  previous  section,  was  one  for  which 
A-^  =  A£  =  -  A  .  This  choice  leaves  three  parameters  to  be  deter¬ 
mined.  They  are  c,  =  A/6  .  and  the  two  constants  a  and  b 

in  the  expression 


A  =  A/\/a  +  b  •  Re^ 

The  method  of  searching  for  values  for  these  parameters  was  as 
follows.  The  equations  for  a  free  jet  were  programmed  so  as  to 
solve  the  system  of  equations  for  a  free  jet  developing  in  the 
axial  direction.  At  an  arbitrary  initial  station  in  the  axial 
direction,  a  mean  velocity  profile  and  profile's  of  the  pertinent 
second-order  correlations  were  arbitrarily  assumed.  For  a  given 
choice  of  model  parameters  (in  this  case,  a  ,  b  ,  and  c.  *  A/r  r  , 
where  r  c  is  the  radius  for  which  u  is  one-half  the  0 

centerline  value),  the  free  jet  equations  were  solved  for  the 
development  of  the  jet  downstream  of  the  initial  distributions. 

In  all  cases,  essentially  self-similar  solutions  were  obtained 
far  downstream  of  the  start  of  the  calculation.  If  a  set  of  para¬ 
meters  could  be  found  so  that  the  resulting  self-similar  flow 
agreed  with  experimental  measurements  with  respect  to  the  rate  of 
spread,  as  well  as  with  respect  to  mean  velocity  and  correlation 
distributions,  it  would  then  be  assumed  that  a  reasonable  tui'bulenee 
model  had  been  achieved. 

Actually,  such  calculations  were  carried  out  for  both  free 
Jets  and  two-dimensional  free  shear  layers.  With  the  single  A 
model,  it  was  found  that  no  combination  of  parameters  a  ,  b  ,  and 
c.  could  produce  an  adequate  description  of  either  a  free  jet  or 
arree  shear  layer.  In  general,  it  was  found  that  if  the  parameters 
were  adjusted  so  as  to  give  an  adequate  rate  of  spread  of  the  mean 
profile  (i.e.,  if  the  level  of  the  turbulent  shear  correlation 
was  large  enough)  the  spread  of  the  correlations  by  diffu¬ 

sion  was  always  too  large.  This  general  result  is  illustrated  in 
Figure  2  where  it  is  seen  that,  if  the  general  level  of  the  shear 
correlation  (u*w’>  were  to  match  the  experimental  data  of 
Wygnanaki  and  Fiedler  (Ref.  11)  in  the  region  cf  maximum  shear, 
it  is  clear  that  far  too  long  a  tail  of  <u‘w*>  at  large  r 
would  result.  This  was  a  very  general  result  for  free  shear  flows 
and  forces  us  to  consider  a  more  complicated  model. 

The  difficulty  that  was  experienced  with  the  constant 
model  was  the  existence  of  too  much  diffusion  relative  to  the  rate 
of  loss  of  correlations,  either  by  dissipation  or  the  tendency 
towards  isotropy.  To  correct  this  difficulty  in  the  studies 
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reported  here,  the  diffusion  lengths  A^  and  A^  were  made 

smaller  than  A^  .  An  idea  of  the  effect  of  reducing  the  diffusion 
lengths  relative  to  the  isotropy  length  can  be  seen  from  Figure  3. 

Here  the  rms  value  of  the  longitudinal  velocity  fluctuation  w’  , 
that  has  been  calculated  for  several  choices  of  model  parameters, 
is  plotted  versus  radius  in  a  self-similar  free  jet.  Note  that  as 
the  diffusion  lengths  A2  and  A^  (which  are  c2  times  A^) 

are  reduced,  the  amount  of  diffusion  is  obviously  reduced  and  the 
levels  of  turbulence  on  the  jet  centerline  are  appreciably  increased. 

The  effect  of  the  choice  of  the  scale  of  the  isotropy  length 
A  can  be  seen  from  Figure  ^ .  The  distribution  of  longitudinal 
turbulence  intensity  is  shown  as  a  function  of  radius  for  two  choices 
of  A^  relative  to  the  local  value  of  r  ^  .  It  is  seen  that  the 

levels  are  much  lower  for  the  smaller  A  than  for  the  larger  value. 
This  is  what  one  might  expect  because  ofthe  increased  dissipation, 
as  well  as  the  increased  loss  of  shear  correlation  by  the  tendency 
towards  isotropy  when  the  scale  A^  and,  hence,  X  is  made  smaller. 

The  effect  of  neglecting  pressure  diffusion  can  be  seen  in 
Figure  5;  the  longitudinal  velocity  fluctuations  in  a  free  jet  are 
shown  as  a  function  of  radial  position  for  a  given  choice  of  model 
parameters  a  ,  b  ,  ,  and  c^  for  two  choices  of  .  One 

choice  is  a  e2  and  the  second  is  «  0,  i.e.,  neglect  of 

pressure  diffusion.  It  is  seen  that  for  this  e no lee  of  the  other 
parameters,  the  effect  of  neglecting  pressure  diffusion  is  not  large. 

Having  given  some  idea  of  how  some  of  the  various  parameters 
entering  the  model  for  turbulent  shear  layers  affect  the  solutions, 
we  must  now  discuss  the  selection  of  an  actual  set  of  parameters. 

If  one  considers  only  a  single  type  of  shear  flow  that  one  wishes 
to  model,  say,  the  free  Jet,  it  is  possible  to  choose  a  whole  spectrum 
of  models  which  will  give  a  good  description  of  the  mean  spread 
of  the  free  jet  and  the  distribution  of,  say,  the  longitudinal 
turbulent  velocity  field.  To  illustrate  this  point,  we  may  refer 
to  Figure  6.  Here  we  see  that  two  profiles  of  longitudinal  velocity 
fluctuation  can  be  obtained  with  radically  different  choices  of  b 
and  A.  .  It  is  observed  that  if  one  chooses  small  b  one  must 
also  enooae  a  small  value  of  A.  relative  to  a  characteristic  scale 
of  the  Jet.  What  then  is  the  basic  difference  between  these  two 
solutions?  It  is  this.  For  the  solution  with  small  b  and  small 
A,  ,  the  balance  of  the  production  of  turbulence  is  mere  by  dissipa¬ 
tion  and  less  by  diffusion  than  for  the  other  ease.  Also,  for  the 
case  of  small  b  and  small  A,  ,  the  solutions  arc  more  isotropic 
on  the  Jet  centerline  than  forAthe  other  case. 

The  choice  between  the  two  models  esMMted  sn  Figure  6  must- 
be  made  on  the  basis  of  the  degree  of  diffusion  and  the  degree  of 
isotropy  desired  in  the  calculated  result.  This  is  a  difficult 
decision  to  make,  for  existing  experimental  data  do  not  agree  as 
to  how  isotropic  free  jets  are  <*w  their  se uteri *r  s,  as  will  be 
seen  presently.  There  is  another  way  that  one  can  decide  between 
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Fig.  3.  Behavior  of  solutions  for  a  seif- similar  free  jet 
when  the  parameter  *  A^  /ti^  8  A^/Aj 
is  varied 


Ig.  4  Effect  of  variation  of  the  isotropy  scale  A 
on  characteristics  of  a  self-similar  free  jet 


Fig.  5,  The  effect  of  neglecting  pressure  diffusion 
when  calculating  a  self-similar  free  jet 


Fig.  6.  Two  choices  of  model  parameters  that  yield 
almost  the  same  distributions  of  <w'w'> 
for  a  self- similar  free  jet 
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two  different  models.  If  one  uses  the  same  model  to  compute  two 
different  turbulent  flows  having  essentially  different  geometries, 
the  model  which  gives  the  best  results  for  both  flows  is,  since 
we  are  seeking  an  invariant  model,  the  one  to  choose. 


As  mentioned  previously,  we  have  computed  self-similar 
solutions  for  a  free  shear  layer  as  well  as  for  an  axially  symmetric 
free  jet.  Actually  a  search  for  model  parameters  for  each  type 
of  flow  was  carried  out.  As  a  result  of  these  studies,  it  was 
determined  that,  insofar  as  the  parameter  studies  have  proceeded 
at  this  point,  the  following  model  for  free  turbulent  shear  flows 
gave  the  best  results: 


a  =  2.5 


b  =  0.125 
=  0.10 
c^  =  0.10 


(26) 


Also,  the  value 


Ax/6 


char 


0.50 


(27) 


was  found  best  for  both  flows,  although  it  was  not  part  oflthe 
plan  to  have  a  common  value  of  cn  .  As  mentioned  above,  for  the 
free  jet. 


char 


(28) 


The  characteristic  length  for  the  free  shear  layer  was  taken  as 


6char  ~  z.25  “  z.75 


(29) 


which  is  the  distance  normal  to  the  flow  in  the  shear  layer  from 
the  point  where  the  velocity  is  one-quarter  the  external  driving 
velocity  to  where  it  is  three-quarters  this  velocity. 

In  Figures  7  through  14,  we  show  comparisons  with  experi¬ 
mental  data  of  the  velocity  correlation  profiles  computed  for  both 
a  free  jet  and  a  free  shear  layer,  using  the  model  parameters  given 
above.  The  experimental  results  are  taken  from  the  work  of 
Wygnanski  and  Fiedler  (Refs.  11  and  12),  Gibson  (Ref.  13),  and 
Donaldson,  Snedeker,  and  Margolis  (Ref.  14). 

Figures  7  and  8  show  the  longitudinal  fluctuations  in  a 
free  jet  and  free  shear  layer,  respectively.  The  agreement  between 
model  calculations  and  experiment  is  good  in  both  cases.  For  the 
free  jet  in  Figure  7,  it  would,  perhaps,  have  been  desirable  to 
have  a  little  more  diffusion  (larger  A,  and  larger  b)  in  the 
model  in  an  attempt  to  reduce  the  oversnoot  in  <Tw,w,'>  near  the 
centerline  of  the  jet. 


Figures  9  and  10  show  distributions  of  normal  fluctuations 
in  both  the  free  jet  and  the  free  shear  layer.  Here  we  note  the 
agreement  with  experimental  data  is  not  so  good.  There  appears 
to  be  a  little  too  much  diffusion  for  these  cases.  Also,  note 
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Fig.  7.  Comparison  of  experimental  results  and  model 
predictions  for  the  longitudinal  velocity  correl¬ 
ations  in  a  free  jet 

<21 


Comparison  of  experimental  results  and  model  predictions  for  the 
longitudinal  velocity  correlations  in  a  free  shear  layer 
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Fig.  9.  Comparison  of  experimental  results  and  model 
predictions  for  the  radial  velocity  fluctuations 
in  a  free  jet 


Comparison  of  experimental  results  and  model  predictions  for  the 
normal  velocity  fluctuations  in  a  free  shear  layer 


Fig.  il.  Comparison  of  experimental  results  and  model 
predictions  for  the  sidewise  fluctuations  in  a 
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Rg.  13.  Comparison  of  experimental  results  and  model 
predictions  for  the  shear  correlation  in  a 


z25~z75 

Fig.  14.  Comparison  of  experimental  results  and  model  computations 
for  the  shear  correlations  in  a  free  shear  layer 


the  very  large  discrepancy  between  measured  normal  fluctuations 
on  the  centerline  as  reported  in  three  separate  experiments.  The 
data  of  Gibson  show  the  components  of  turbulent  velocity  to  be 
essentially  isotropic  on  the  jet  centerline,  while  those  of 
Wygnanski  and  Fiedler  and  Donaldson,  Snedeker  and  Margolis  do  not. 
From  the  results  shown  in  Figure  9,  it  would  appear  that  if  one 
were  to  desire  more  isotropy,  one  would  wish  to  choose  a  smaller 
value  of  A.  and,  hence,  a  smaller  value  of  b  .  This  is  opposite 
to  the  conclusion  drawn  from  Figure  7. 

Figures  11  and  12  show  the  sidewise  components  of  turbulence 
for  the  free  jet  and  free  shear  layer,  respectively.  The  agreement 
between  experiment  and  computed  results  is  better  for  the  free 
jet  than  for  the  free  shear  layer.  The  reason  for  this  behavior 
is  not  known. 

In  Figures  13  and  14,  we  show  the  shear  correlations  for 
the  free  jet  and  the  free  shear  layer.  The  agreement  in  both  cases 
is  fair.  It  should  be  noted  that  the  experimental  values  of  shear 
correlation  from  Ref.  11  have  been  shown  as  reported  (the  open 
circles)  in  Figure  14  and  also  as  corrected  hy  us  (the  solid  symbols) 
so  as  to  agree  with  the  measured  rate  of  spread  of  the  free  shear 
layer.  A  comparison  of  the  measured  shear  and  that  inferred  from 
the  mean  velocity  profile  was  reported  by  Wygnanski  and  Fiedler 
but  apparently  their  computations  contained  an  error.  Also  shown 
in  Figure  14  is  the  lev:l  of  shear  that  may  be  inferred  from  the 
mean  spread  of  the  free  shear  layer  studied  by  Tollmien  (Ref.  19) 
and  frandtl  (Ref,  It?)  many  years  ago.  It  is  seen  from  the  results 
presented  in  Figures  13  and  14  that  the  model  gives  a  fairly  good 
representation  ©f  the  shear  in  both  the  free  jet  and  the  free  shear 
layer. 


A  easeful  study  of  Figures  7  through  14  shews  that  it  really 
is  necessary  to  study  further  the  prsMem  of  choice  of  model  para¬ 
meters.  However,  before  this  is  done,  it  appears  desirable  to 
have  at  hand  experimental  data  which  one  can  rely  on  to  be  truly 
representative  of  the  basic  flow  which  is  being  calculated.  It 
is  difficult  to  choose  a  more  sophisticated  model  until  the  question 
of  the  degree  of  isotropy  on  the  centerline  of  a  free  jet  is  settled. 
In  addition,  one  should,  at  this  point,  determine  if  the  model 
just  found  for  free  shear  layers  can  be  used  for  a  model  of  the 
ou^er  regions  of  a  boundary  layer  and  give  reasonable  results. 

8-sfore  turning  to  the  problem  of  the  turbulent  boundary  layer. 
It  will  be  instructive  to  find  a  relationship  between  the  values 
of  A,  used  in  the  free  shear  layer  and  the  free  jet  calculations 
and  the  general  magnitude  of  the  Integral  scales  measured  for  such 
flows.  In  the  computations  that  have  been  made,  it  has  been  assumed 
that  A  i§  constant  across  a  free  jet  or  a  free  shear  layer  at 
any  given  longitudinal  position  and,  in  magnitude,  proportional 
to  the  local  scale  of  the  mean  flow.  It  is  well  known  that  the 
integral  scales  of  such  flows  are,  in  general,  proportional  to  the 
local  mean  scales  but  the  actual  value  of  the  integral  scale  varies 
across  the  layer. 


In  Table  1,  we  present  the  values  of  integral  scale  within 
a  free  jet,  as  reported  by  Wygnanski  and  Fiedler.  The  integral 
scale  tabulated  is  the  longitudinal  integral  scale 


L  = 


<w '  w  ’>  (z 


<w 1 


(J21)WI 


(2o)>  ^(Zo  ■  ZT  ) 


(30) 


for  the  free  jet. 


Table  1 


Integral  Scales  in  a  Free  Jet  after  Wygnanski  and  Fiedler  (Ref.  11) 


Radian  Position  Dimensionless  Scale  Scale  Ratio 


r/x 

0 

.05 

.10 

.15 

,20 

Also  presented  in  Table  1  is  the 
A  ,  to  the  local  integral  scale 

this  ratio  for  the  free  jet  Is 
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Q.448 
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0.595 
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0.59 
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Table  * 

integral  .Scaled  in  a  Free  Shear  Layer  after  ifygnan^k*  and  Fiedler 

%*er.  UJ 

Location  in  Jet  Dissent  ion  less  Scales  Scale  Satie 


tr/X 

L/fs.25"2.75) 

V*“ 

inner  Region 

0 . 093 

0.«*5 

0.59 

Center 

0.103 

0.883 

0.8? 

Outer  Region 

0.14? 

2  .2? 

0.39 
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A  typical  value  of  A,/L  for  a  free  shear  layer  appears  to  be 
approximately 

Ax/L  =  0.55  (33) 

As  mentioned  previously,  there  is  not  much  point  in  going 
further  with  studies  of  the  present  model  until  it  has  been  applied 
to  a  boundary  layer.  Note  should  be  made  at  this  point,  however, 
of  other  methods  of  calculating  turbulent  shear  flows  -  methods 
that  are  similar  to  the  methods  being  discussed  here.  As  mentioned 
before,  the  idea  behind  the  method  is  not  new.  It  follows  a  trend 
suggested  by  Prandtl  and  Wieghardt  (Ref.  1)  and  follows  closely 
the  general  line  of  approach  taken  by  Rotta  (Ref.  6).  Since  these 
two  pioneering  papers,  many  more  or  less  similar  studies  have  been 
undertaken.  Typical  of  many  such  studies  are  those  of  Glushko 
(Ref.  18),  Bradshaw,  Ferriss  &  Atwell  (Ref.  19),  Harlow  &  Romero 
(Ref. 20),  Gawain  &  Pritchett  (Ref.  21),  and  Beckwith  &  Bushnell 
(Ref.  22).  The  more  sophisticated  of  these  methods  do  not  assume 
a  local  scale,  as  we  have,  but  carry  along  with  the  computational 
scheme  an  equation  for  the  required  scale.  Such  an  equation  was 
derived  by  Rotta  (Ref.  6)  from  the  equation  for  the  two-point 
correlation  tensor  <u?(A)uj(B)>  .  There  is  no  question  but  that, 

in  the  future,  the  method  presented  here  should  be  enhanced  by 
coupling  the  present  set  of  model  equations  to  an  equation  for  the 
integral  scale.  To  date,  however,  we  have  avoided  making  this 
connection  in  order  to  study  the  character  of  the  model  and  its 
dependence  on  the  scale  A-,  without  this  dependence  becoming  inex¬ 
tricably  mixed  with  the  additional  modeling  that  must  be  done  in 
the  equation  that  is  used  to  compute  a  scale. 

Before  proceeding  further,  it  must  be  demonstrated  that, 
if  the  present  model  is  applied  to  a  boundary  layer,  useful  results 
will  be  obtained  for  the  same  choice  of  model  parameters  that  has 
been  made  for  free  turbulent  shear  flows. 


V.  APPLICATION  OP  MODEL  TO  BOUNDARY  LAYERS 


If  the  model  of  turbulent  shear  flows  is  to  be  applied  to 
a  boundary  layer,  the  parameters  ,  c^  ,  a  and  b  are  known. 

But,  since  the  characteristic  length  in  a  boundary  layer  is  arbitrary 
(as  it  is  in  the  free  jet  and  the  free  shear  layer),  we  are  at 
liberty  to  choose  c,  ,  i.e.,  the  ratio  between  A,  and  the  charac¬ 
teristic  length  (which,  in  this  case,  we  take  equal  to  6  ^  ,  the 

thickness  cf  the  layer  in  which  the  velocity  reaches  99%  of  its 
free  stream  value). 


As  discussed  in  a  previous  section  of  this  paper,  one  other 
parameter  enters  the  problem,  namely,  a  ,  the  coefficient  appearing 
in  Eq.(21).  We  have,  then, 


=  aVa  z 

for  0  z  c-^6  ^^/(av/a)  ,  and 

A1  =  c16.99 


(34) 


(35) 


for  z  >  c16  ^/(on/a) 


With  only  these  two  parameters  a  and  c,  to  determine,  the  search 
is  not  difficult.  The  model  that  has  been  found  is  the  following: 


a  =  2.5 
b  =  0.125 

c2  =  ® 
c3  =  ® 

=  Ax/6>99  =  0.15 
a  =  0 .7//~a  =  0 . Lh3 


(36) 


The  ability  of  this  model  of  a  turbulent  shear  layer  to 
predict  the  known  mean  properties  of  turbulent  boundary  layers 
is  shown  in  Figures  15  through  17.  In  Figure  15,  we  show  the  skin 
friction  developed  by  our  model  as  it  proceeds  from  a  disturbed 
laminar  layer  to  a  fully  turbulent  layer.  Also  shown  are  the  lami¬ 
nar  skin  friction  law  and  the  turbulent  law  proposed  by  Coles 
(Ref.  23)  which  is  a  good  fit  to  experimental  data.  It  is  no  great 
surprise  that  the  general  levels  of  skin  friction  we  computed  agree 
well  with  experimental  findings  inasmuch  as  the  values  of  a  and 
c,  were  chosen  to  get  those  levels  correct.  Of  more  importance 
is  the  nearly  exact  following  of  the  trend  of  skin  friction  with 
Reynolds  number  by  the  model  computations. 
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Fig.  15.  Computed  variation  of  coefficient  of  friction  with  Reynolds 
number 
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Fig.  16.  Computed  velocity  profiles  for  three  Reynolds  numbers. 

Computation  started  a*  Rex  =20,000.  The  reference  velocity 
uT  is  defined  by  rwa||=/»u|  . 
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Fig.  17.  Computed  velocity  defects  for  three  Reynolds  numbers. 

Computation  started  at  Rex=20,000.  The  reference  velocity 
uT  is  defined  by  rWOjj sPur  * 


Figure  16  shows  a  comparison  of  the  computed  mean  velocity 
profiles  developed  by  the  model  In  the  vicinity  of  the  wall  and 
the  well-known  law  of  the  wall  as  proposed  by  Coles  (Ref.  23). 

It  may  be  seen  that  the  law  of  the  wall  is  not  quite  achieved  by 
the  present  selection  of  model  parameters.  The  results,  however, 
are  sufficiently  accurate  to  be  encouraging. 

In  Figure  17  we  compare  the  experimentally  determined  velocity 
defect  law  proposed  by  Coles  (Ref.  23)  with  the  results  of  our 
model  computations.  It  is  seen  that,  once  the  turbulent  boundary 
layer  is  well  established,  the  computational  model  gives  a  fairly 
good  representation  of  the  outer  regions  of  the  turbulent  boundary 
layer. 


With  these  results  in  hand,  we  must  now  consider  the  relation¬ 
ship  of  the  computational  scales  used  to  the  longitudinal  integral 
scales  that  are  found  in  the  outer  regions  of  turbulent  boundary 
layers.  For  this  purpose,  we  may  use  the  measurements  of  Grant 
(Fief.  24).  When  the  experimental  correlations  reported  by  Grant 
for  y/6  =  0.66  *  in  a  turbulent  boundary  layer  are  integrated  to 

give  the°longitudinal  integral  scale,  one  obtains  L/6o  »  0.3  . 

Since  for  our  calculations,  6  /6  =»  0.33,  we  find  that 

L/6  ^  a»  0.25  •  Since  the  computational  scale  used  was  A^/6  ^  = 

Q.15j  we  find  that 

VL  -  °-“  iw 


This  is  a  moat  welcome  result  sine©  it  shows  that,  for  all  the 
turbulent  flows  we  have  investigated,  the  ratio  of  the  proper  compu¬ 
tational  scale  to  the  longitudinal  integral  scale  is  approximately 
the  same. 


•Grant  defined  6^  as  that  height  in  the  bmmdary  layer  where  the 
velocity  defect  iSw  was  equal  to  the  friction  velocity. 
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VI.  CALCULATION  OF  ATMOSPHERIC  TURBULENCE 


We  may  now  turn  to  our  original  problem,  namely,  that  of 
computing  atmospheric  turbulence.  In  Ref.  2,  uhe  equations  for 
the  generation  of  atmospheric  turbulence  were  given  and  computa¬ 
tions  were  made  for  some  simple  shear  layers,  using  the  originally 
determined  model  parameters  [Eqs.(24)].  The  mode]  equations  for 
the  calculation  of  atmospheric  turbulence  using  the  latest 
modeling  are  given  in  Appendix  D. 

In  what  follows,  we  will  turn  our  attention  to  a  calculation 
of  a  case  for  which  atmospheric  turbulence  has  been  observed  and 
measured  by  an  instrumented  aircraft  and  for  which  data  were 
taken  on  the  mean  state  of  the  atmosphere  nearby,  at  a  time 
reasonably  close  to  the  time  of  the  observation. 

In  searching  the  records  available  from  the  LO-LOCAT  Phase 
III  atmospheric  turbulence  study,  it  was  decided  that  an  appropri¬ 
ate  LO-LOCAT  run  on  which  to  make  our  study  was  leg  5  of  Test  b3 
(Category  *122224).  In  Figure  IS  we  reproduce  from  Ref.  25  the 
normalized  atmospheric  spectrum  determined  from  this  particu¬ 

lar  leg  of  Test  43.  The  spectrum  is  complete  enough  so  that  both 
the  turbulent  intensity  and  a  longitudinal  integral  scale  can  be 
estimated.  We  shall  treat  this  spectrum  as  a  partial  spectrum, 

i.e.,  as  if  only  the  4*  ( Q )  —  -  power  portion  of  this  spectrum 

was  available.  Using  only  this  information  and  the  information 
available  from  the  state  of  the  atmosphere ,  we  will  attempt  to 

estimate  the  values  of  and  L  by  calculation  of 

the  turbulence  that  can  be  produced  by  the  atmosphere  for  various 
values  of  the  scale  L  . 

The  log  for  this  flight  leg  on  20  September  1968  shows 
(Ref.  2b): 


Time 

Airspeed  (frue) 

Altitude 

Ambient  temperature 
9 round  temperature 
Wind  a.. eed 

Wind  direction  (true) 
at  150  feet 

The  state  of  the  winds  at  the 
19. 


Mid -morning 
6 1 4  ,Y  feet/sec. 

715. 7  feet 

9  2.7*F 

33.5  feet/sec 

light  altitude  is  shown  in  Figure 


Aircraft-  615  ft/sec 
heading  348°  true 


Wind*  33.5  ft/sec 
from  169°  true 


Rg.  19-  Wind  and  aircraft  velocities  at  flight 
altitude,  Test  43,  Leg  5  (Ref.  26) 


Meteorological  data  gathered  at  0900  hours  local  standard  time 
on  the  same  day  at  a  point  about  15  miles  from  the  flight  leg  are 
available  v  ief .  27).  These  data  have  been  plotted  in  Figure  20; 
shown  are  mean  velocity  and  mean  temperature  measurements  deduced 
from  the  data,  with  lines  faired  through  the  points  to  define  an 
atmospheric  mean  velocity  and  mean  temperature  profile.  Also 
shown  are  the  ground  temperature  and  ambient  fat  750  feet)  temper¬ 
ature  and  wind  velocity  as  recorded  in  the  flight  test  log.  The 
considerable  discrepancies  are  to  be  expected,  given  the  distance 
and  time  intervening. 

In  Figure  21  the  temperature  profile  has  been  replotted  to 
show  the  departure  of  the  local  temperature  from  an  adiabatic 
temperature  profile.  It  will  be  seen  that  the  temperature  profile 
is  such  that  it  is  unstable  from  the  ground  to  a  height  of  some 
900  feet;  thereafter  the  temperature  profile  is  stable  in  the 
region  of  interest. 

Using  these  profiles,  the  equations  for  the  development  of 
atmospheric  turbulence  presented  in  Appendix  D  can  be  solved  for 
the  turbulence  that  will  be  generated.  The  only  free  parameter 
in  the  model  is  the  computational  scale  .  From  the  discussion 
presented  in  previous  sections,  we  will  assume,  for  these  calcula¬ 
tions,  that  this  scale  behaves  as  follows. 

Near  the  ground,  the  scale  is  given  by  the  formula 

A  *  a  SEz  =3  0.72 

It  will  behave  in  this  manner  until  it  reaches  some  arbitrarily 
assumed  upper  limit  A£  .  Above  this  altitude,  A1  is  assumed  to 

be  constant  with  altitude  and  equal  to  Aif  . 

Although  we  know  that  the  integral  scale  and,  hence,  the 
computational  scale  generally  increase  with  altitude,  it  is 

believed  that,  since  we  are  interested  here  in  a  computation  of 
turbulence  in  only  the  lower  levels  of  the  atmosphere,  the  assump¬ 
tion  of  A  constant  above  the  level  of  large  ground  interaction 
will  be  adequate. 

Computations  were  carried  out  for  a  range  of  values  of  Af 
equal  to  60,  100,  200,  300,  and  400  feet.  In  Figures  22  and 
23,  the  profiles  of  \u'u')  ,<v‘v,s>  ,  <w'wf>  ,  and  the  shear 
correlation  <u'w*>  are  plotted  as  a  function  of  altitude  for 
values  of  At*  of  100  and  200  feet,  respectively.  For  completeness, 

we  also  show  in  Figures  24  and  25  plots  of  the  profiles  of  the 
temperature  flux  correlations  <u*T'>  and  s^w'T^  as  well  as  values 

2 

of  the  temperature  variance  <T'  )  . 
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Fig.  20.  Reported  wind  and  temperature  values  from 
Ref.  271  along  with  test  leg  values  (Ref.  26) 
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Fig, 24.  Computed  profiles  of  temperature  correlations 
for  A*  «  IOO  ft 
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Rg.  25.  Computed  profiles  of  temperature  correlations 
for  A*  =  200  ft 
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If,  from  results  such  as  those  plotted  in  Figures  22  and 

23,  one  determines  the  values  of  o  =  <(u 'u ')> ,  0  =  <(  v  ' v ’>  , 

1  /p  ^  v 

and  c?w  =  ‘(w'w'y  x  at  the  flight  altitude  (750  feet),  a  plot  of 

these  values  as  functions  of  the  outer  integral  scale  A^  can  be 
constructed.  Such  a  plot  is  shown  in  Figure  26.  The  very 
pronounced  effect  of  the  scale  A*  on  the  magnitude  of  the  turbu¬ 
lence  that  can  be  developed  is  exhibited. 

We  must  now  compare  these  results  with  the  spectral  results 
already  discussed.  In  comparing  the  three  values  of  A-i/L. 
determined  previously,  namely,  0.69  for  the  free  jet,  0.55  for  the 
free  shear  layer,  and  0.60  for  the  boundary  layer,  we  select  the 
relationship 


A1  =  A*  =  0.60L 


as  representative. 

In  Figure  27  we  show  c?w  again  as  a  function  of  A|  as  in 
Figure  26  and  as  a  function  of  L  .  Also  shown  is  a  curve  of  the 

I/O 

relationship  between  o  and  L  [o  ~  L  ,  according  to  Eq . 
(5)]  obtained  from  the  Spectrum  shown  in  Figure  18,  treated  as 
a  partial  spectrum.  We  note  that  the  two  curves  intersect  at 
a  =  3.0  ft/sec  and  L  =  290  feet.  This  may  be  compared  with  the 
results  obtained  from  the  complete  spectrum  (Ref.  25),  namely, 
ct  =  2.81!  ft/sec  and  L  =  273  feet. 
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The  remarkable  agreement  between  these  results  must  be 
considered  fortuitous  in  view  of  the  many  uncertainties  involved 
in  constructing  wind  and  temperature  profiles  for  input  to  our 
program  from  the  few  measurements  points  available  and,  more 
particularly,  in  assuming  these  values  apply  some  15  miles  away 
at  a  different  time.  The  discrepancies  shown  in  Figure  20 
between  the  two  sets  of  measurements  give  an  indication  of  the 
validity  of  our  results. 

Although  the  agreement  we  have  just  found  is  better  than 
expected  and,  certainly,  additional  calculations  of  the  type  we 
have  just  presented  must  be  carried  out,  we  do  believe  we  have 
demonstrated  a  method  of  computing  atmospheric  turbulence  that 
has  sufficient  merit  to  warrant  further  investigation  and 
development.  The  other  conclusion  that  one  can  draw  from  the 
calculation  that  we  have  pi'esented  is  that,  in  the  future,  it 
would  be  most  desirable  to  expend  every  effort  to  obtain  measure¬ 
ments  of  the  mean  state  of  the  atmosphere  as  close  in  time  and 
space  to  any  actual  flight  measurements  of  turbulence  as  is 
possible . 


L ,  ft 

Fig.  27  Calculated  turbulent  intensity,  ,  as  a  function  of 
compared  with  the  relationship  derived  from  the 
spectrum  according  to  Eg.  (5) 


VII.  DISCUSSION  AND  CONCLUSIONS 


We  have  in  this  report  presented  a  method  for  calculating 
turbulent  shear  flows  that  depends  on  the  closure  of  the 
equations  of  motion  at  the  second-order  correlations  of  fluctuating 
quantities.  This  closure  is  accomplished  by  modeling  the  appropriate 
terms  in  the  equations  for  the  second-order  correlations  in  terms 
of  simple  invariant  expressions  involving  the  second-order 
correlations  themselves  and  several  scale  parameters.  By  comparison 
of  computations  with  experiment,  the  relationships  between  the 
various  scales  were  determined  and  the  relationship  of  these 
scales  to  the  integral  scale  in  the  turbulent  flow  that  was  being 
calculated  was  found.  It  was  found  that  the  basic  scale  used 

in  the  computation  was  approximately  equal  to  three-fifths  of 
the  longitudinal  integral  scale  of  the  turbulence  under  consideration. 

Although  it  should  be  pointed  out  that  more  computations  must 
be  made  before  a  final  choice  of  model  and  model  parameters  is 
made,  it  is  clear  that  some  of  the  inconsistencies  in  the  presently 
available  experimental  data  on  free  turbulent  shear  layers  should 
be  cleared  up  first. 

Once  a  turbulence  model  had  been  established  and  a  relation¬ 
ship  found  between  the  computational  scale  used  and  the  integral 
scale  of  the  turbulent  motion  under  consideration,  the  method  was 
applied  to  the  problem  of  obtaining,  from  partiai  atmospheric 
spectra,  an  estimate  of  the  total  turbulent  intensities  and  the 
integral  scales.  8y  applying  the  method  to  a  case  of  low  altitude 
turbulence  where  the  spectrum  was  complete  enough  to  determine 
both  the  intensity  and  scale  of  turbulence  but  treating  the 
spectrum  as  if  it  were  incomplete,  it  was  found  that  the  magnitude 
of  the  turbulent  intensity  and  scale  could  be  evaluated  if  the  mean 
wind  and  temperature  profiles  were  known. 


Although  the  method  appears  to  work  well,  further  studies  of 
this  type  should  be  performed.  At  the  present  time,  it  is  felt 
that  the  technique  that  has  been  developed  is  of  sufficient  interest 
to  warrant  support  for  these  additional  studies.  An  important  part 
of  these  studies  would  be  an  application  of  the  method  to  high 
altitude  turbulence  where  the  magniUide  of  the  integral  scale  is 
not  known. 


In  connection  with  such  studies  of  high  altitude  turbulence, 
it  is  strongly  recommended  that  all  possible  effort  be  expended  to 
obtain  mean  atmospheric  wind  and  temperature  data  as  close  in 
space  and  time  to  the  actual  measurement  of  turbulence  as  is 
possible . 
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APPENDIX  A.  CONSTANT  DENSITY  MODEL  EQUATIONS  FOR  STEADY  FLOW 


Continuity : 


u^,  =  0 
>  J 


(A.l) 


Momentum: 


'  -  P  5,1  *  '  (eJkSi.k).J  -<u')'ui>,J 


Stress  tensor: 


sJ<uiuk>,j =  -<uJ'uk>  5u 

♦  gJ*  [A2q(<uiuj>  k  +  <u]u’>jl  *  J 

*  [a3o  <“r“i>,i],k  *  |V  i],i 


-  ^  -««!)♦  “tJ£<“i“k>,4X  -  2 


APPENDIX  B. 


MODEL  EQUATIONS  FOR  A  BOUNDARY  LAYER 


We  use  a  Cartesian  coordinate  system  (x,y,z)  with  the  free 
stream  velocity  in  the  x-direction  and  with  a  as  the  coordinate 
normal  to  the  wall.  The  velocity  components  are  denoted  by 
(u,v,w).  With  the  usual  assumptions  for  a  constant  pressure 
boundary  layer,  Eqs.  (A.l)  through  (A. 3)  become 


Bu  .  dw  „ 

33T'  SF  3  0 


(B.l) 

(B.2) 

(B.3> 


-  Bu  .  -  Bu  _  ..  B^u  B<u‘w'> 

^  ^  - 

u  ~  +  w  »  -2<u»w'>  | “  +  ( M 

Bx  B2  32  02  \  *  32  / 

-  — -  I) 


(B.3) 
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u  sr — -  +  3  1 z  )  *  i^(<v  v  >  -  j) 
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(B.M 


-  ^»'w')  .  r.  B<V'w*>  _  B  r 
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(3Aj  ♦  **j)q  -H<yw>- 

J  X 

<*.s> 

u  ♦  s  . .  <«.„.>  I  ♦  F(2A2  ♦  A3)q 
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APPENDIX  C.  MODEL  EQUATIONS  FOR  A  FREE  JET 

We  use  the  cylindrical  coordinate  system  (r,0,z)  with  the 
velocity  components  denoted  by  (u,v,w),  respectively.  The 
assumptions  for  a  constant  pressure  boundary  layer  are  again 
applied  to  Eqs.  (A.l)  through  (A.3)>  along  with  the  assumption 
of  axial  symmetry,  giving 


|i+ 1  +  |*  =  o 

or  r  bz 


a^+-^w/afw  +  l^w  \  _  _  <u'w’> 


dr 


dz 


L^r 


2  r  Sr 


dr 
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♦£( 
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-  d<v'v'>  ,  -  d<v'v'>  _  d  /  a  „  d<VfV’>  \  A  3A2q  S<v'v’> 
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-  0<W,W,>  ,  -  3<w'w' 

u  +  w  —  ■  - - - 

dr  dz 
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APPENDIX  D.  MODEL  EQUATIONS  FOR  CLEAR  AIR  TURBULENCE 


As  in  Ref.  2,  we  simplify  the  problem  somewhat  by  assuming 
the  mean  flow  is  in  the  x  direction  only  and  that  all  the 
dependent  variables  are  functions  only  of  time  t  and  altitude 
z  .  The  convective  operators  in  the  steady  state  equations 
(U/cix.j )  are  replaced  by  S/St  . 


In  considering  atmospheric  motions,  we  include  buoyancy  terms 
and  temperature  terms  that  induce  them.  We  take  the  actual  mean 
temperature  to  be  the  sum  of  an_adiabatic  temperature  T0(z)  and 
a  departure  therefrom  denoted  T(z,t).  The  density  corresponding 
to  the  adiabatic  temperature  variation  is  denoted  by  p0(z)  and 
we  write  n 

P0 

The  equations  before  modeling  are  given  in  Ref.  2. 
current  modeling,  they  become 


With  the 


aS  32u  a(P0<u'"'>) 
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In  these  equations,  a*  is  the  second  coefficient  of  viscosity 
and  g  is  the  acceleration  of  gravity. 
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